Improving approximate matrix factorizations for implicit time integration in air pollution modelling
نویسندگان
چکیده
For a long time operator splitting was the only computationally feasible way of implicit time integration in large scale Air Pollution Models. A recently proposed attractive alternative is Rosenbrock schemes combined with Approximate Matrix Factorization (AMF). With AMF, linear systems arising in implicit time stepping are solved approximately in such a way that the overall computational costs per time step are not higher than those of splitting methods. We propose and discuss two new variants of AMF. The rst one is aimed at yet a further reduction of costs as compared with conventional AMF. The second variant of AMF provides in certain circumstances a better approximation to the inverse of the linear system matrix than standard AMF and requires the same computational work. 2000 Mathematics Subject Classi cation: 65M12, 65F30, 65M20, 65M06 1998 ACM Computing Classi cation System: G.1.8, G.1.3, J.2
منابع مشابه
Numerical Time Integration for Air Pollution Models
Due to the large number of chemical species and the three space dimensions, o -the-shelf sti ODE integrators are not feasible for the numerical time integration of sti systems of advection-diffusion-reaction equations @c @t +r (uc) = r (Kr c) +R (c) ; c = c(x; t); c 2 IRm; x 2 IR3 from the eld of air pollution modelling. This has led to the use of special time integration techniques. This paper...
متن کاملA New Approximate Matrix Factorization for Implicit Time Integration in Air Pollution Modeling
Implicit time stepping typically requires solution of one or several linear systems with a matrix I − τJ per time step where J is the Jacobian matrix. If solution of these systems is expensive, replacing I−τJ with its approximate matrix factorization (AMF) (I − τR)(I − τV ), R+V = J , often leads to a good compromise between stability and accuracy of the time integration on the one hand and its...
متن کاملSolving vertical transport and chemistry in air pollution models
For the time integration of stii transport-chemistry problems from air pollution modelling, standard ODE solvers are not feasible due to the large number of species and the 3D nature. The popular alternative, standard operator splitting, introduces artiicial transients for short-lived species. This complicates the chemistry solution, easily causing large errors for such species. In the framewor...
متن کاملOn the Efficient Treatment of Vertical Mixing and Chemistry in Air Pollution Modelling
Reaction-Advection-Diffusion problems occurring in Air Pollution Models have very stiff reaction terms and moderately stiff vertical mixing (diffusion and cloud transport) terms which both should be integrated implicitly in time. Standard implicit time stepping is by far too expensive for this type of problems whereas widely used splitting techniques may lead to unacceptably large errors for fa...
متن کاملILU and IUL factorizations obtained from forward and backward factored approximate inverse algorithms
In this paper, an efficient dropping criterion has been used to compute the IUL factorization obtained from Backward Factored APproximate INVerse (BFAPINV) and ILU factorization obtained from Forward Factored APproximate INVerse (FFAPINV) algorithms. We use different drop tolerance parameters to compute the preconditioners. To study the effect of such a dropping on the quality of the ILU ...
متن کامل